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Abstract 

We consider a one-parameter family of non-evolutionary partial differential equa- 
tions which includes the integrable Camassa-Holm equation and a new integrable equa- 
tion first isolated by Degasperis and Procesi. A Lagrangian and Hamiltonian formu- 
lation is presented for the whole family of equations, and we discuss how this fits into 
a bi-Hamiltonian framework in the integrable cases. The Hamiltonian dynamics of 
peakons and some other special finite-dimensional reductions are also described. 



£j 1 Introduction 

> 



In this note we consider the following one-parameter family of partial differential equations 
OO : (PDEs): 

o ■ 

o ■ 

On 
O 

CN • version of the integrable Camassa-Holm equation, which is 

o ■ 



Ut - Uxxt + (& + = hUxUxx + UUn^xx (1-1) 

(the parameter h is constant). In the particular case b = 2 this becomes the dispersionless 



Ut ~ Uxxt + 3UU^ — 2UxUxx + UUxxx- (1-2) 

With the inclusion of linear dispersion terms the equation (|1.2| ) takes the form 

XXX ^xxt '^^x^xx ~l~ ^^xxx (■^■'^) 

(co,7 and a are constant parameters), and in this form it was derived as an approximation 
to the incompressible Euler equations, and found to be completely integrable with a Lax 
pair and associated bi-Hamiltonian structure 0. 

Degasperis and Procesi 0] applied the method of asymptotic integrability to test a many- 
parameter family of equations generalizing ( |1.3| ), and found that only three equations passed 
the test up to third order in the asymptotic expansion, namely KdV, Camassa-Holm and 
one new equation. After rescaling and applying a Galilean transformation, the new equation 
may be written in dispersionless form as 

Ut - Uxxt + 4:UUx = SUxUxx + UUxxx- (1-4) 
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Henceforth we refer to ( |1.4| ) as the Degasperis-Procesi equation, and observe that it is the 
6 = 3 case of the equation In a recent article @] we proved the integrabihty of the 

Degasperis-Procesi equation by providing a Lax pair, and derived two infinite sequences 
of conservation laws. The key to our construction was to use a reciprocal transformation 
connecting the equation ( |1.4| ) with a negative flow in the Kaup-Kupershmidt hierarchy. 
Both the Camassa-Holm equation (|L^) and the Degasperis-Procesi equation ( |1.4| ) display 



the weak Painleve property |T^, with algebraic branching in their solutions. Thus when 
Gilson and Pickering applied the standard Painleve tests to a class of equations including 
the family (|1 . 1| ) , these two integrable equations were excluded because of the branching, 
and they found that no equations passed the tests. In a forthcoming work ||^ we show that 
after applying a reciprocal transformation to the family it is possible to use Painleve 
analysis effectively, and the requirement of only pole singularities immediately isolates the 
two integrable cases b = 2,3. Because the equations (|1 . 1|) are non-evolutionary, the standard 



symmetry approach of Shabat et al |^ does not apply. However, recently Mikhailov and 



Novikov developed a powerful extension of the symmetry classification method [T^, and 
applying this to the equations (|TTT|) they found that only the cases 6 = 2,3 could possess 
infinitely many commuting symmetries, and so only these two cases are integrable. 

One of the novel properties of the Camassa-Holm equation is that it admits exact solutions 
in terms of a superposition of an arbitrary number of peakons, or peaked solitons. The A^- 
peakon solution is just 



N 



«(x,t) = ^p,(t)e-l--'^^WI, (1.5) 
i=i 

where the positions qj{t) of the peaks and their momenta Pj{t) satisfy an associated dynam- 
ical system taking the canonical Hamiltonian form 

dh . dh 

with 

j,k=i 

It is shown in [|l], 0| that this is a completely integrable finite-dimensional system, with an 
N X N Lax pair. The integrabihty of the Camassa-Holm peakon dynamics further received 
an r-matrix interpretation in |jTl] , 



One of our observations in was that all of the equations in the family ( |L1|) admit 
A^-peakon solutions in the form ( |L5| ), but the dynamical system for qj, pj only takes the 
canonical Hamiltonian form (|1.6| ) in the particular case 6 = 2. In the following we construct 
a Lagrangian (for 6 7^ 0) and associated Hamiltonian structure, depending on the parameter 
6, for each equation in the family (|1.1|) . We then reduce this onto the finite-dimensional sub- 
manifold of peakon solutions for each N to obtain a new non-canonical Poisson structure and 
associated Hamiltonian form for the peakon dynamics. It turns out the two-body dynamics 
(A^ = 2) is Liouville integrable for any value of 6, but for arbitrary we expect integrabihty 
only for 6 = 2,3 since these are the reductions of integrable PDEs. An N x N matrix Lax 
matrix pair for the 6 = 3 peakon dynamics was presented in 0], but at that stage we were 
still ignorant of the correct Poisson structure. 



2 Lagrangian form and Hamiltonian operator 

In the following we consider solutions of ( p..l|) on the whole a;-axis vanishing at infinity, and 
all integrals will be assumed to be between ±00, although (with suitable modification) most 
of the results hold for solutions with constant background or on a periodic domain. Some of 
the operations below only make sense for sufficiently smooth real functions u, but later we 
will want to apply some formulae to weak solutions such as the peakons, in which case we 
will find that certain manipulations are valid also for distributions. 
It is convenient to introduce the quantity 

m = u — Uxx, 

which is just the Helmholtz operator acting on u, so that the family of equations ( |1 . 1| ) may 
be rewritten in the form 

rrit + urrix + buxfn = 0. (2-8) 
For any value of b this is a conservation law, i.e. 

mt + (^^^^-^^{u^ - ul) + urn^ =0, 

while for all 6 7^ there are at least three conserved quantities, namely 

J mdx, j m^/^dx, j m-^'^ (^-^ + 1^ dx. (2.9) 



For the integrable cases 6 = 2 [Q, |^, and 6 = 3 |Q there are infinitely many conserved 
quantities. 



For the second quantity appearing in ( ^.9] ) above, the associated conservation law is 
conveniently written as 

Pt + ipu)x = 0, m = -p\ (2.10) 

With this further rewriting of the equation ( |2.8D , or equivalently ( |1.1| ), for 6 7^ we introduce 
a potential rj such that 

Vx =P, Vt = -pu. 

Then in terms of this potential we find that (apart from the particular cases b = 0,1) the 
equation (LI) may be derived from the following action: 

S = J j Cdxdt = j j Q^[(logr/,.)xx + l] + ^^dxrft. (2.11) 

When 6 = 1 a similar Lagrangian formulation is valid, but the Lagrange density must be 
modified to ^ 

^ = 2^tVx^ [(^ogrix)xx + l]+Vx ^ogr]x. 

Starting from the Lagrangian we can apply a Legendre transformation in the usual way. 
The conjugate momentum is 

C = Tr = ^ [(logr/x)zx + 1] , 
dr]t 2r]x 

and (for b ^ 1) the Hamiltonian is 

H= {(rit- C)dx = --^ / r]^^dx = -^ / mdx. (2.12) 



Having applied the Legendre transformation we then find that for any b the PDE ( |2.8| ) can 
be written in Hamiltonian form as 

mt = B—, (2.13) 



where the operator 



B = —{bmdx + rnx){dx — d^) {bmdx + {b — l)m. 



(2.14) 



is skew-symmetric and satisfies the Jacobi identity (for a proof see (P). In the case 6 = 1 
the Legendre transformation gives a different result: J mdx is a. Casimir for the operator B, 
and instead the Hamiltonian may be taken as 



H 



m log m dx. 



The formula ( p.l3|) is even valid in the case 6 = 0, when the Lagrangian formulation breaks 
down. 



For the integrable case 6 = 2, the Camassa-Holm equation (|1.2| ), there are two local 
Hamiltonian structures [ffl, 13] , given by 



Bo 



-d.{i-dl), B, 



and the compatible operator ( p. 14 ) is the first nonlocal Hamiltonian structure obtained by 
applying the recursion operator R = BiBq^ to Bi, i.e. 

B2 = BiBq^Bi = B\b=2- 

In the other integrable case 6 = 3, namely the Degasperis-Procesi equation (|1.4|) , there is 



only one local Hamiltonian structure, and the operator ( |2.14|) gives the second Hamiltonian 
structure, viz 

Bo = -9,(1 - dl)iA - dl), fii = fi|,=3. (2.15) 

This compatible bi-Hamiltonian pair first appeared in [^]. The nonlocal part of (p.l4|) may 
be defined precisely via 



1 r°° 

{dx-dl)-'f{x) = -j G{x~y)f{y) 



dy, 



G{x) = sgn{x){l 



-\x\ 



)• 



(2.16) 



The Jacobi identity for (|2.14]) reduces to a single functional equation satisfied by the kernel 

G m. 



3 Peakon dynamics and special solutions 



The peakon solutions of ( |1.1| ) take the form ( |1.5| ) for any 6, in which case the dependent 



variable m is a sum of Dirac delta functions, 

N 



m{x,t) = 2j2 Pj it)Si^ - Qj (t) ) • (3.17) 



Introducing the even kernel (Green's function of the Helmholtz operator) 



g{x) = e-^^\, G{x)= I g{s) ds 



X 







the dynamical system for the peakons may be conveniently written as 



(3.18) 

^ = -{b- 1) J2k=iPjPk9'i(lj - Ik)- 

Clearly this takes the canonical Hamiltonian form ( |1.6|) only for 6 = 2. 

In order to obtain the correct Poisson structure on the reduced peakon phase space, for 
any 6 7^ 1 we can use the Poisson bracket defined by (|2.14] ), which gives (up to rescaling) 



{m{x),m{y)} = _ (G{x - y)m^{x)m^{y) + hG'{x - y)m{x)m^{y) 

— hG'{x — y)m{y)m^{x) — b'^G"{x — y)ra{x)ra{y)^ . (3.19) 

Substituting the peakon expression (|3.17| ) into (p.l9|) , we are able to calculate the Poisson 
bracket on the reduced (2A^-dimensional) phase space as 

{Pj^Pk] = -{b - l)G"{qj - qk)PjPk, {qj,Pk} = PkG'iqj - Qk), 

{q„qk} = l/{b~l)G{q,-qk). (3.20) 

The Jacobi identity for this non-canonical Poisson bracket follows from the same functional 
equation for G that arises for the operator (|2.14|) (see [^). The Hamiltonian for the dynamical 
system ( |3.18| ) is found simply by substituting ( |3.17] ) into the PDE Hamiltonian ( |2.12| ), so 
that after scaling we have 

N 

h = {b-l)H/2 = J2Pr 
i=i 

Having obtained the reduction of the Poisson structure for the PDE, we see that the 
peakon dynamical system ( p. 18 ) takes the Hamiltonian form 



dt dt 
Now for any b we also introduce the second quantity 

N 



^ = \^ PjPk(^ - (1 - 9{(lj - (lk)f ^ 
j,k=i 



For 6 = 2 this just coincides with the canonical Hamiltonian ( |1.7| ) for Camassa-Holm 
peakons, k = h, and is clearly conserved. For the Degasperis-Procesi case 6 = 3, with the 
bracket ( ^.20 ) k Poisson commutes with h for any A^; it is the second conserved quantity 



obtained from trace L^, where L is the Lax matrix presented in However, for any other 
value of 6 it seems that k is conserved only in the two-body case N = 2. Thus we see that 
for generic values of the parameter 6 only the two peakon dynamics is Liouville integrable. 
The two-body dynamics and associated peakon phase shifts are described in more detail in 

The integrability of the two-body dynamics seems to induce remarkable stability proper- 
ties for the peakons even in the non-integrable regime. Numerical studies of the PDE family 
( p..l| ) on a periodic domain 0] show that starting from a Gaussian initial profile for u, for any 
6 > 1 a train of peakons emerge and preserve their shape and speed after interaction, under- 
going only a phase shift. The peakons are unstable for 6 < 1. For the range — 1 < 6 < 1 the 



numerical results indicate that "ramps-and-cliffs" profiles, like those appearing as solutions 
to the Burgers equation, seem to be stable. We do not have an exact analytic expression for 
such profiles. However, the "ramp" part behaves like the exact similarity solution 



while for the special case 6 = we find an exact solution in terms of a superposition of N 
"cliffs," 

N 

u{x, t) = c + ajG{x — ct — Qj{t)), 
i=i 

with c and the aj being constants. The Qj satisfy the Hamiltonian dynamical system 

dQ- . . 

= akG{Qj - Qk) = {Qj, h}, ^ = "j'^J' 
Mi i=i 

where (after scaling) the Poisson bracket reduced from ( p.l9| ) is 

{Qj, Qk} = G{Qj — Qk)- 

The "cliff" dynamics is integrable for N = 2 and = 3. 

In the remaining range 6 < — 1, the numerical results of show that Gaussian initial data 
splits into a superposition of pulses which are asymptotically stationary. For that parameter 
range there is a unique vanishing stationary solution given by the exact sech-shaped profile 

, 2/(b+l) 

u{x) = i cosh[(6 + l){x — Xq)/2] ' 



centred around the arbitrary point xq. In future work we intend to address the stability 
properties of the peakons and other exact solutions. 
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